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1. INTRODUCTION 

Since the rediscovery of von Zeipel's method by D. Brouwer (1959) 
and its successful application to the problem of artificial satellite s , many 
other problems have been solved by that same method, thus proving 
its great applicability. It is the purpose of these notes to present the 
general equations of Von Zeipel's method and discuss briefly their 
applicability. 

The reduction of the order of a differential canonical system can, 
in theory, be performed by obtaining, one way or another, integrals of 
the system. One of them is the Hamiltonian itself when it is time 
independent. Actually, this integral of the system (physically its 
"energy"), can describe completely the geometry of the solutions in a 
phase space of 2n dimensions where 2n is the order of the system. 
When this order is 2, then the solution is completely specified and the 
use of the Hamiltonian reduces it to a first order differential equation 
which can be integrated by quadrature. The introduction of integrals 
in a system of n degrees of freedom (2nth order), reduces it to one of 
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n -p degrees of freedom which can be integrated immediately when 
n - p - 1 (where, of course, p cannot be greater than n). 

A few comments can be made with respect to the more famous 
methods of reduction to show their eventual connection with von Zeipel's 
method* 


2. FROM HAMILTON TO von ZEIPEL 

In the discussion that follows only methods that have been used in 
connection with differential systems describing the motion of a physical 
system are considered. The presentation does not necessarily follow a 
chronological order. 


Consider then a system of n_ degrees of freedom given by 2n first 
order canonical equation 


BH 



(j = 1,2, . ■ ■ ,n) (1) 

BH 



where the Hamiltonian H = H (x , ..., x n , y 2 , ...y n )is presumed to be 
time independent. If this is not the case, the introduction of time as a 
new canonical coordinate x n+1 (the associated momentum being -H ) 
always reduces the latter to the former case. The degree of freedom 
will however increase by one. 
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A canonical transformation of the variables (x,y) to new variables 
(x', y') will be, in this exposition, equivalent to the problem of finding 
a generating function S = S (x f , y, t ) such that 



(i = 1, 2, . - - , n) 


x. 

j 



( 2 ) 


It is easily seen that this is a sufficient condition to satisy the Jacobi- 
Poincare relation 


(Xj dy . - x'. dy'p = dW . (3) 

j = i 

which is valid whether or not S is an explicit function of time. The 
Hamiltonian of the new system will be equal to that of the old one inas- 
much as one is obtained from the other by introducing the transforma- 
tion of variables expressed by Equations (2) when 3S/3t = 0. 

a. HAMILTON- JACOBI — The method introduced by Hamilton and 
Jacobi consists in obtaining a canonical transformation such that the 
new Hamiltonian is identically zero. In such a case, the new variables 
are all constants. 

b. LIND5TEDT ’S METHOD — Linds tedt's method is a particular 
application of the Hamilton- Jacobi method for cases where the 
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Hamiltonian is expanded in terms of small parameters. In this partic- 
ular case the solution gives the coordinates as linear functions of time 
and the momenta as constants (usually called action angle variables). 
The comparison with the Hamilton- Jacobi method is purely formal 
since the method devised by Lindstedt is quite original. Actually, the 
real difference between von Zeipel's and this method is that Lindstedt 
does not make use of a generating function. 

c. WHITTAKER r S METHOD (solution by series). This method 
obtains n integrals of the system by reducing the Hamiltonian to a 
function of the products p j = x. y, ( j = 1,2,... ,n). In this case, since 

3H m 

~ - Wt Xj 

bh m 

y i “ Bx. Bp. y j’ 

J *J 

it follows that 

x. y. + x. y . = 0 or p. = const (j = 1, 2, ... ,n). 

d. DELAUNAY'S METHOD — This method, as Lindstedt's, can be 
applied only when the Hamiltonian consists of a "zero order" part (the 
corresponding system having a known solution) and a "disturbing 
function" that has a small numerical factor. The basic approach of the 
von ZeipeLs method is the same as that of Delaunay’s method; how- 
ever, the latter one makes no use of a generating function and breaks 
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the disturbing function into parts which are treated separately. The 
Hamiltonian must be constructed after the transformation is performed 
for each particular part. 

A few more techniques could be mentioned but one deserves more 
attention than all the others. The concept of adiabatic invariants in 
Quantum Mechanics is quite analogous to the concept of "mean variables 1 ' 
in von Zeipel's method, or to a certain extent to what Whittaker calls 
Adelphic Integrals. 

3. THE von ZEIPEL'S METHOD (1916) 

It has been quite common, after Delaunay, to use the negative of 
the Hamiltonian. Thus, if F=-H and if 4. (j = 1 ,2,...,n) and L. (j = 

1, 2, ..., n) are the coordinates and momenta respectively, then 


dF 

dL, 


(j =1,2 n) 


L. 


) 


iL 

) 


(4) 


Suppose 


F =F(^, L; e) (5) 

where e is a "small parameter" and Z and Lindicate the sets , ..., ^n) and 

(Lj L n ). A canonical transformation involving the parameter e will 

be given by a generating function 
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S = set, L*; €) 


such that 


t _ 3S 

i " 3-t. 


(j =1,2 n) 


( 6 ) 


C =2^. 
j 


where (-£*, L* ) are the new coordinates and momenta. If F* is the 
negative of the new Hamiltonian, then we assume 

F* (l\ L*; e) = L*; e), L (%* , L*; e); e) (7) 


or, from Equations (6), 


F* 




') 


In a more restrictive sense it is assumed that the series 


h= n 


e k F k (l, L) 


( 8 ) 


( 9 ) 


represents the negative of the Hamiltonian to the required degree of 
precision and converges to F(l, L; e) as F rom this point F is 

written as F without danger of confusion. Furthermore, it is assumed 
that 
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S(Z, L*; e) 


F k (l, L) 

F*(l\ L*; e) 

are developable in Taylor's series in the neighborhood of e = 0, so 
that the series 


CD 

k = 0 



are convergent for sufficiently small 6. 


By the conservation property 


L 

n 



(ii) 


where it is important to note that BS/B't contains e through Equa- 
tions (6). Equating the coefficients of like powers in e in both sides of 
Equation (11) gives a system of partial differential equations in Sand 
F*. The next step is obtaining this system. 

4. DIFFERENTIAL EQUATIONS OF THE VON ZEIPEL'S METHOD 
The derivative of F k with respect to e at the point e = 0 is 
obtained as follows. 
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Consider 


dF k _ yi ^F k dLj y-> BF^ ^ / BS\ 

de / i 3L de L . — i de ^3-t. J 


Using Equation (10) it follows that 

/dFA y, 3F k d J yi (-t.L-)' l 

\ de j /L 3L, de ) ^ f 

= Y T j ei-> (I?i) . 

4 - 4 3L > 

Let us now compute 

de " 1 " 1 V 3 L i/ 


( 12 ) 


where 


F, = F, 



Applying Leibniz' formula, this becomes 


l m- 1 


\ min (m~ 1 , j~l ) 

¥ ■ L M 

1 1 v=0 


:J-1 


de 


m- 1 


i-l\ d v e* -1 d’ 


m-l-v 


de" de 


m-l — v 



For e = 0 the only possible choice is j <m. Then 



(J 
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Now using Equation (12) 


fdm F k 

k de m 


n cd 




€-° i=1 j=l 



It is now desirable to rewrite the above equation as 


fd m F, 


^de m / e=0 


n oo 


■LL>> 

i 1=1 i 1=1 


(m - 1)! / BS i i 




d i i /-d'F. 


de m_i l \ 3L n 


6 = 0 


(13) 


Equation (13) is now applied to find 


d m ~M l“ x k 


d€ m "Ji V i 



e=o 


The result is 


d — i. 


de m- ii \^ij 


J-i J-i . On-jj - 1)! 
‘ 2 = 1 i2 = 


BS . 


12 


i— i L-i 2 (m-j -j )! 

; ° i,= i j o=l 1 1 x 


2 / e = 0 


d m ~ii i 2 

\1 

de -n-il-i 2 

K 3L * j_ 


e=0 


The process is repeated up to the point where 


m -J, -J 2 - 


j N - 0. 


(14) 


so that 


!-i 2 iN 


B H F, 


j m — J 1 — J2 — - • - — JN \BL. BL. . . . BL 

de J1 Si 1N \ n ‘2 inJ 


6 = 0 


b n f, 


BL. . . . BL . , 

l l 'N/e^O 
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Substituting these successive derivatives into Equation (13), it follows 


that 


n co 


M m F. 


de m / £ =o i 1= i j 1 = i 


(m-1)! it, x 


1 (m - jj)'- l 


*1/ e=0 


n 00 


y-i v — i . (m-jj-1)! / 3S 

Z_I Z_I (tn - jj - j 2 )! 

i 2 =l j 2 =i v 


J 2 


2/6=0 


n co 


>3 =1 J3 =1 


. C"> - j , - j, -1)! / ' 3S i 3 '\ x 

(m-j.-jj-j,)! ^,3 A, 0 


n ® 


■zz: •• 

X N =1 JN =I 


(m - j i - j 2 - • • • Jn-i ^ ^ 


JN 


(m - j , - j 2 j N )! 


l N / 


e=0 


B N F, 


BL. BL. ■■•BL. . 

J i >2 *N/ £ =o 


The numerical factors are reduced to 


h 


m! 


N m(m - j j) (m- jj - j 2 ) "•(m- jj - j 2 j N _j) 


£ m! C(m; j x , j 2 , ' * * , j N ) 


and the second summation does not run in general up to infinity but to 
a limit given by condition (14). Thus the above relation becomes 


© -zii (£)««*> '-'lift) ( 


B n F, 


BL. ...BL, , 

*1 1 N/ e=0 

(15) 
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where ^ stan ds for summation over all possible positive integers _j p 

(m) N 

whose sum is tn (according to Equation (14)), and the first product J | 

’ " p=i 

refers to the summation signs • There are n of these integers. 

1 p = 1 

Equation (15) will be valid even for m = 0 with the definitions 


d Fu 


de 


- F, 


O k 


and 


C(0; -) = 1. 

From Equations (10) and (11) it follows that 


UJ CD CD / 1 m T~' \ 

k = 0 k-0 m=0 ' 06 1 6 = 0 


00 CD 


k-0 m=0 


em+k ^d m F k \ 

m! 

lde m ) 

jse 

results 


CD v 


L L^r 


v /d m F 


e = 0 v - 0 m= 0 


£ = 0 


The substitution of these results into Equation (14) leads to 

’•ttzfitty (=%-) •• 

^=0 m=0 (m) p=l \i p =l/ P =1 \ 1 P/e = 0 \ 1 1 I N/ / f = o 


(16) 


In a complete similar way, if 

00 

f, =£ £kF ;(I .' ■ v Y 


k=0 
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then 


'■LLLU (L 

v-Q m=0 (m) p=l yip -1 

(17) 

It is important to note that in Equation (16), e - 0 is equivalent to 
L r = 3S Q /d't r ( r = 1, 2, and in Equation (17) e- 0 is equivalent 

to i* = 3S Q /3L* (r = 1, 2 n), according to Equation (10). The 

equality of factors of the same power of e in Equations (16) and (17) 
gives the partial differential equations for the von Zeipel's method 


C(m; j r 


-J 


-JV/3S. 

■>n*r 



3^* 

>N / ,-n 



C(m; jj, 




+ 


3S. 

J p 




= 0 


for v = 0,1,2,... . 


(18) 


For instance, Equation (18) gives: 


v = 0 


f „ e, 


!!£ 

34 




(19) 


v = 1 



1 
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'1 


05 OJ 

^ o W 


a; a> 

-H- 00 

QJ | 05 

in 


05 

f-t QJ 

— KJ 


a5 a 5 

^ vi 

XT O 


Q J 05 

<±o in 


05 | Q 5 

r ^ 


051 QJ 

C/3 

H- o 


sM 


05 05 

r c/3 

►•• *|m 

a> | a> 


a; a; 

f C/3 
“• ■» O 


05 Qj 

r* c/5 

* o 




05 | QJ 

r.|j" 

oj | o5 

^3 

M- #1 t-* * 


cy 05 

qj\ a > 
c^> I C/3 

o5 

jr sd 

05 

r ° 


05 05 

<r^ C/3 

o 


05 % 


sM- 


05 05 

05 I 05 

r tj 


5^ 05 

xj\ 


QJ QJ 

r c /3 

* o 


sM- 


05 05 

r* c/j 

M- # H- 


05 Q5 

in 

M* H* 

05 05 

in 


05 Q5 

*H- W 

05 

M P 05 
05 w 

jr <s* 

oJ 

r 


05 05 

-Hk, W 

051 05 

r ^ 


05 | 05 

^o| C/3 


05 05 

r in 

~+ \ ~ 

05 | 05 

r w 


<r^ Kj 
M ' * 

05 o * 
<cN> 


05 I 05 

r L w 


05 05 

r w 

* o 




05 05 

in 


05 05 

r *3 

— o 


05 05 

cn 

— o 


05 05 

<^o *^3 


05 05 

r « 


*0 

II 

05 I 

05 

c^l 

c/3 

■o | 

o 


+ 




where use has been made of the coefficients 
C (1; 1) = 1 

C (2; 2) = 1 C (2; 1, 1) =y 

C (3; 3) = 1 C (3; 2, 1) =-| C (3; 1, 2) =y C (3; 1, 1, 1) =-|- . 

5 . elimination of variables 

Since the solution of the system is known where F is reduced to F Q , 
the problem is to eliminate variables which are not present in F Q . 
Suppose a canonical transformation is found in such a way that_p of the 
n coordinates (p < n ) have been eliminated from the Hamiltonian, that 
is 
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(2 3 ) 


f = r ' k*. *;■ l ;. l ; l ;) • 


The equations of motion then yield 


L* = C k (const.) (k = 1, 2, . . . , p) . 


If these constants are replaced in F*, then 


F* 



V. C,. c 2 . 



and the problem is reduced to one of n-p degrees of freedom. 


a. If p = n , the problem is completely solved, since 



(k = 1, 2, n) 


and 


K = 


(Cj, 


C n ) t + V k (0). (k = 1, 2, .... n) 


b. If p = n-1, the problem is integrable by quadrature. In fact 


L* = C k (k = 1, 2, .... n-1) 



n 


Since 


F (Cj, . . . , C n l ; ^ n , L n ) - C - const. , 


(24) 
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then 


L 


* 

n 


K (c. <V 



and therefore 


V = X 

n n 


(C. C, c n _, : r) 


and 


t - 


t o ” 


.r 


dC 


k*(0) ( C > C X C n— 1 ’ O 


The coordinate ■i* becomes a known function of time as well as L* 

n n 

Therefore, the equations 


K 


ar 

^L k * 


\ (c,. C 2 , . 


- c n =i; l : (t). ^(t>) 


(k = 1, 2, 


n-1 ) 


can be integrated by quadrature. 

The von Zeipel’s method consists in the elimination of some of 
the coordinates (angular variables) and the reduction of the problem 
to case (b) and possibly (a). The adaptability of this method is based 
on a set of hypotheses which are listed below in Roman numerals. 

I) The new and old corresponding variables differ by a quantity at 
least of the first order, i.e. 
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0 (€) 


(i = 1, 2 n) 

L* - L. = 0(e). 

This automatically fixes S Q to correspond to the identity trans- 
formation since for e = 0, the above conditions give 

V. = l. 

1 1 

(i = 1, 2, . . . , n) 

L* = L. . 

1 1 

Therefore 


s.=2y l :- (25) 

irl 

If expression (25) is substituted into Equations (19), (20), (21) 
and (22) , then 

v = 0 


F 0 (4, L*) 

V = 1 

F, {l, L*) 
= F* (4, L*) 


f* c-e, L*) 



(26) 


(27) 
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V = 2 





and similarly for Equation (22). 

It is seen that S k (L*, 1) and F* (L*, -C) are unknown functions. In 
order to perform a particular solution toward the elimination of certain 
angular variables in F* we impose conditions (which are usually suitable 
in Celestial Mechanics) on the functions S k and F* . They are 

II) F* (L*. l) does not depend on <£. (i =1, 2, .... p < n) for any 
k > 0. 

Ill) S k (L* , l) only depends on the 4 . ( i = 1 . 2, .... n) through 

trigonometric functions, for any k > 0. This avoids "secular 
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perturbations” in the momenta L . , or in other words differences 


L 

j 



B (S - S Q ) 
) 


are periodic functions of the (k = 1, 2, n). 


The application of these conditions, together with the obvious fact 
that F 0 does not depend on angular variables 'C. ( i = 1 , 2, p < n ) 
which are to be eliminated, yields the relations 


+0 (< 


P+1’ 


l 


K) = K (*,„■ 


K) <«> 


F 


1 5 


= K 


F 


ip 




( 30 ) 


F 


2s 


+ P 


2s 



+ P 


♦ 

2 s 


F 2p + P 2p + 

i = 1 





( 31 ) 


and so forth. The functions F, and F, , F 0 and F 0 , P 0 and , P* 

Is Ip 7 2s 2p 2s 2p 2s 

and P* p are the portions of Fj, F 2 , P 2 and P* which are respectively 
independent of and dependent on the { i = 1, 2, p) , and where 
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p 


2 




3 2 Fp 

BL* 3L* 

i i 


P 


* 

2 “ 





a 2 K 

M. M. ' 

i J 


( 32 ) 


In the usual problems of Celestial Mechanics F Q does not depend on any 
angular variable so that the Equations (30), (31), (32) and the corre- 
sponding equations for higher order are much simplified. Thus, the 
additional hypotheses will be considered. 

IV) F 0 and thus F* depend only on the momenta L* 

V) The angular variables <1 (i = 1, 2, ..., m) corresponding to 
momenta L x (i = 1, 2, . m ) that are present in F Q have been 
eliminated to the k th order. 


The next problem is the possibility of elimination of angular variables 
whose conjugate momenta are not present in F Q . At this stage the 
Hamiltonian of the system is 


f - = f; (l; l;) + f; (*„, 




+ 



.... t n . l;. 



where 

L* = C = const (j - 1, 2, m) , 

j j 7 


( 33 ) 
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and the old and new variables are related by 


3S. 3S, 

t T * _ £ £ 

i " i " dt + U. 

J J 


3S k 

-U. 


(j = 1, 2 n) 


(34) 


as, 

V. - i. = — - 

J J BL* 


bs 2 

oh* 

1 


!!i. 

3L* 

J 


Assume a new canonical transformation from the variables 




m + 1 9 


..., V , L 


m+ 1 


L**^ and let 


L*) to the variables (£**., 

n / V m+ 1 



L 


* * 

m+ 1 ’ 


s' - s* (r v, l;;, l;-) 


(35) 


be its generating function. Then, since L* = Cb = const ( j = 1,2,..., 

m). 


75* / T ** T **\ * / T ** T **\ 

0 ( L i - L m ) = F o ( L i L J = 


- const 


(36) 


L* k * = L* = C k = const (k = 1, 2, . . . , m) 


f ; (^,1 K - c r c 2 c „. c, l v) 

•r k k -- l ;;. o ■ 


(37) 


The last equation implies that the elimination of further variables is 
possible if and only if F* does not depend on them. For in this case 
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i (K 


m + p+ 1 ’ 


V, c, 


'm+1’ 


■ C) 


= F 


^ l ;:, o 


and 


bs: bf: 


v — 1 os> or. 

F! (-T r ; L*\ L*‘) + ) - - 

2 V m+I’ n m+1 n / / J 3^* BL* 


i = m+ 1 


V-. BS* 3 FV 

= fv (r r ; l , , .... v*) + > — - 1 

2 \ m+p+1’ n’ m + 1 n ) BL** B^* 

i =m+ p+ 1 i * 


(38) 


which defines S* . It is important to note that in such a case S* will be 
defined by an equation involving; 2nd order terms; these terms are 
therefore necessary to obtain first order ’’perturbations . M This fact is 
exactly what happens in Brouwer’s theory on artificial satellites (1959), 
where 


a) The elimination of g*is possible because F* is independent of 
this variable. 


b) The development for long period perturbations ( those of argument 
g*) needs the evaluation of 2nd order terms. 

This type of reasoning can be carried on up to any order in exactly 
the same way. It may then happen that the elimination of a certain 
angular variable by obtaining Sj requires the evaluation of terms of the 
kth order. 
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However if F* depends on the angular variables to be eliminated 
the problem cannot be solved unless it happens that the remaining 
system has one degree of freedom. For example, this is the case of 
the perturbations on the motion of an artificial satellite by the moon. 

6. SMALL DIVISORS 

The case of critical inclination for the theory of artificial satellites 
of an oblate planet for which P 2 is the dominant zonal harmonic and 
J 4 / is a well known example of the problem of small divisors. 

Here, only a particular aspect of the question is dealt with. Consider 
the solution of Equation (30) in the usual case where F* does not depend 
on the 't . The characteristic associated system is 


d'Cj 

d^ 2 

<u 

P _ 

dSj 

dF o 

' dF o " " 


F lP 

bl; 

bl; 

3L* 

P 



(39) 


Should one of the partials dF 0 /dL* happen to be zero, the general solu- 
tion would certainly be discontinuous since a n small divisor" is present. 
However this divisor is not exactly zero because the quantity 3F Q /3L* 
is evaluated to first order only. 

In the case of critical inclination it is necessary to take 

S = S 0 + e 1/2 S 1/2 + c S 1 + e 3/2 S 3/2 + ... - 

However, in doing so the separation of "long periodic' 1 and "secular" 
perturbations is lost. The integration leads, in most cases, to elliptic 
integrals (Hori, I960). 
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The question of small divisors usually arises whenever the problem 
presents cases of libration as particular solutions. 

Another case to be mentioned is the resonance for an artificial 
satellite whose period is commensurable with the period of rotation of 
the Earth when tesseral harmonics are included. Again, expansion in 
powers of e 1/2 can be used to solve the problem (Morando, 1962). 

Finally it is important to note that singularities in the Equations 
(39) reflect singular points in the hypersurface defined by the Hamil- 
tonian of the system in a phase-space of 2(n-p) dimensions if p variables 
have already been eliminated. 

7. SUMMARY 

The general differential equations of the von Zeipel's method have 
been given to any order. It is hoped that this will avoid tedious Taylor 
expansions if one needs to go to order higher than the second. 

At the same time, the brief discussion on the applicability and a 
few pathological cases of the method, may give some guidance toward 
the solution of new problems. 
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